Local Moment Formation in the Periodic Anderson Model with Superconducting 
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We study local moment formation in the presence of superconducting correlations among the 
/—electrons in the periodic Anderson model. Local moments form if the Coulomb interaction 
U > U cr - We find that U cr is considerably stronger in the presence of superconducting correlations 
than in the non-superconducting system. Our study is done for various values of the /—level energy 
and electronic density. The smallest critical U cr values occur for the case where the number of /— 
electrons per site is equal to one. In the presence of d— wave superconducting correlations we find 
that local moment formation presents a quantum phase transition as function of pressure. This 
quantum phase transition separates a region where local moments and d— wave superconductivity 
coexist from another region characterized by a superconducting ground state with no local moments. 
We discuss the possible relevance of these results to experimental studies of the competition between 
magnetic order and superconductivity in CeCu2Si2- 

PACS numbers: 75.20.Hr, 71.27. +a, 74.70. Tx 
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The superconducting and magnetic properties of 
heavy-fermion materials have attracted much attention 
mainly due to their non-conventional characterBEl All 
these materials have very large specific heat coefficients 
7, indicating very large effective masses, hence the desig- 
nation heavy fermions. The complexity of these systems 
arises from the interplay between Kondo screening of lo- 
cal moments, the antiferromagnetic (RKKY) interaction 
between the moments and superconducting correlations 
between the heavy quasi-particles. On the experimental 
side, these systems exhibit phases in which antiferromag- 
netic ordering of the local moments may coexist with un- 
conventional superconductivity, and/or phases with no 
magnetic ordering. The latter could be either due to 
Kondo screening of the local moments or to a spin liquid 
arrangement of the local moments. 

For systems that exhibit both superconductivity and 
antiferromagnetism, usually [/-based heavy-fermions, the 
ratio between the Neel temperature T/v and the super- 
conducting critical temperature T c is of the order of 
Tn /T c ~ 1 — 10, with coexistence of both types of order 
below T c . The coexistence of both types of order can be 
tuned by external parameters such as external pressure 
or changes in the stoichiometry.EpQ Examples of heavy- 
fermion materials which exhibit antiferromagnetic and 
superconducting order at low temperature are URU2S12 
and Uo.97Tho.o3Bei3. It has recently been found that 
UPd 2 Al 3 (T N = 14.3 K and T c = 2 K) and UNi 2 Al 3 
(Tjv = 4.5 K and T c = 1.2 K) show coexistence of sr 
conductivity and local moment antiferromagnetism.! 
In the Ce-based heavy-fermions however typically mag- 
netism competes with superconductivity. In the proto- 



type heavy-fermion system Cea;Cu2Si2 the competition 
of d— wave superconductivity and magnetic order was 
clearly identified in a small range of x values around 
x ~ 0.99.1 

A description of the normal state properties of the 
heavy-fermion systems has been attempted assuming a 
generalisation of the impurity Anderson model to the 
latticexnl3 In the Anderson lattice the energy of a single 
electron in an /—orbital (e. g. 4/ 1 ) is eo, and the energy 
of two electrons in the same /—orbital (4/ 2 ) is 2eo + U, 
where U is the on-site Coulomb repulsion. The energy of 
the 4/ 2 state is much larger than the energy of the 4/ 1 
state. In the case of interest the empty /—orbital lies 
below the Fermi level ej?, whereas the doubly occupied 
/—orbital has higher energy than ep. Since in many sys- 
tems eo + U 3> £f the simplifying limit U — > 00 is taken 
is some theoretical studies. 

The limit E/p=f p? has been studied using the slave bo- 
son technique. E2rt2l In particular, it has been shown that 
superconducting instabilities arise in the p and d-wave 
channels because of thepoffective (RKKY) interaction be- 
tween the /— electrons .113113 Mean-field studies of super- 
conductivity in the Anderson latti ce , bqt h at finite U and 
at U = 00 have recently been done£3'Ej Nevertheless, the 
coexistence of superconducting correlations and magnetic 
ordering of local moments in heavy-fermion systems has 
not yet, to our knowledge, been studied theoretically. 

In this work we aim to establish the conditions 
for the appearance of local moments in the Anderson 
lattice when superconducting correlations among the 
/—electrons are also present. In order to decide whether 
a /-site behaves as a local moment we use a criterion 
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identical to that introduced by Anderson for the single 
impurity problemli3,r-which we shall now briefly review. 
The procedure in Re£a consists of making a Hartree-Fock 
decoupling of the interaction at the impurity (/ site): 

Uh^nf -> J7(n{)n{ + f/n((n{) - U(n\)(fi[) , (1) 

where , n| are the number operators at the / site 

and their expectation values (nl),(h{) must be found 
self-consistently. Then, the symmetry beaking solutions 
with m = (hi) — (h{) 7^ correspond to the local mo- 
ment regime. Therefore, the local moments arise as the 
symmetry-breaking local minima, with (hi) ^ (^f)) of 
the effective action for the impurity in a Hartree-Fock 
decoupling scheme. Such a mean-field treatment (Q) does 
not account for the dynamics of the local moment which 
arises from its effective interaction with the conduction 
electrons. 

In this work we propose to extend such ideas to the 
lattice of /-sites taking phenomenologically into account 
the presence of superconducting correlations between the 
/-electrons. Because we consider a lattice of / sites, we 
have to look for mean-field solutions with some previ- 
ously chosen spatial arrangement of the moments which 
will be taken as either ferromagnetic or antiferromag- 
netic. Since heavy fermion systems tend to exhibit an- 
tiferromagnetic order, we have chosen this type of bro- 
ken symmetry state. Some comments are also made on 
the ferrromagnetic case. We recall that this mean-field 
treatment leaves out the dynamics of the local moments. 
Such dynamics lies beyond the scope of this work but 
would possibly allow to describe the nature of the mag- 
netic phases of the system including Kondo screening, 
or ordering due to RKKY interaction, or a spin liquid 
arrangement of the local moments, etc. 

Our study of the interplay between superconducting 
correlations and local moment formation may be relevant 
to the understanding of recent studies of CeCu2Si2 sam- 
ples near stoichiometric composition, where it has been 
observed that a ci-wave superconducting phas e exp els a 
magnetic "A phase" under increasing pressureEZHLS. 

We wish to investigate the appearance of local mo- 
ments in the Anderson lattice in the presence of super- 
conducting correlations. Because a microscopic descrip- 
tion of superconductivity in heavy- fermion systems is still 
lacking we opt to treat superconductivity by explicitly 
adding a phenomenological pairing term Hj to the An- 
derson lattice Hamiltonian for spin 1/2 electrons: 

H = H Q + Hu + Hj, (2) 

where 

+ ^E( c Ua-+4^) > ( 3 ) 

£,<T 



H u = U^h lA h lA , (4) 

i 

H ' = j E(4/-m a w + ILc -) > ( 5 ) 

ft 

where e^ is the dispersion of the c— electrons, eo is the 
bare energy of the localized /—states, V is the hy- 
bridization matrix element (assumed k independent), 
U is the on-site Coulomb interaction, and rij jCr = 
fiafi.o- The gap function A(fc) can be written as 
r)(k)A, where rj(k) indicates the pairing symmetry and 
A = jj- Yl%(f-% \ f% j) 7 ^- To simplify the calculations we 
consider a two dimensional system. The results for three 
dimensions are similar. In 2D we have^jj = cos(fe x ) + 

cos (fc^), vf — sin(fci), »7g = cos(fe x ) — cos(k y ), for 
s, p, and d wave, respectively. 

The local moment formation is investigated by apply- 
ing the decoupling (jlj) to the term Hu- We seek mean- 
field solutions such that the electron occupation at the 
/-site i is given by 

(flafi^) = \ n f+ \<rmcos(Ri ■ Q), (6) 

where n / and Ri denote the mean number of / electrons 
and the position of the i lattice site, respectively, and 
a = ±1. This choice corresponds to a spatially periodic 
arrangement of the local moments where Q is the order- 
ing vector. Since we shall perform the calculation on a 
2D square lattice the choice Q = gives ferromagnetic 
order and the choice Q — (tt, ir) gives antifcrromagnetic 
order. 

The mean-field parameters m, rif = n| + n', and A 
are obtained from the minimization of the free energy. 
The chemical potential is determined imposing total par- 
ticle number conservation n = n c + rif . We perform the 
calculations for a conduction band dispersion of form: 
eg = — 2£ J2 i=x cos(fcj) . All the energies and tempera- 
tures in this work are measured in units of the hopping 
integral t, and we have set t = 1. The possible pairing 
symmetries have been studied separately. 

For a given set of parameters n, eo, V, and J one may 
ask what is the minimum value of the local repulsion, 
U cr , for which ground-state solutions with (A, m) =/= 
appear. Figure [j] shows l/U cr for each eo and two elec- 
tronic densities. The local moment regime appears for 
U > U cr . For comparison we also show U cr for a non- 
superconducting system (J = 0). It is seen from Figure 
|l| that the appearance of local moments in a supercon- 
ducting system requires considerably stronger U than in 
a normal system. For each curve the / level occupancy 
rif increases upon increasing eo and the values of eo at 
which U cr is minimum (left panel of Figure |l|) correspond 
to 71/ = 1 . Therefore we also see that smaller U cr occurs 
when / level occupancy is close to one. If rif < 1 then 
U can be increased without destroying superconductivity 
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but if rif > 1 theu-iacreasing U will eventually destroy 
superconductivity.E3'li3 We note that the relative position 
of the superconducting curves in the two panels of Figure 
|l| is different. 




FIG. 1. Relation between the critical U cr and the / level 
bare energy €q, for the parameters J = —1.7, V = 0.4, t = 1, 
and: left panel: n = 2; right panel: n = 1. The curve for 
the normal system corresponds to J = 0. Calculations were 
performed for temperature T = 10~ 3 . 

We have also considered the case of a ferromagnetic 
arrangement of the local moments (Q = 0). Our calcula- 
tions show that, for the same parameters, ferromagnetic 
solutions require much stronger U values (1/U cr < 0.16) 
even in the absence of superconductivity. This seems to 
be an indication that, in a lattice, there is a tendency for 
local moments to form in an antiferromagnetic arrange- 
ment. 

The dependence of the Neel and superconducting tem- 
peratures on pressure has been measured in some heavy- 
fermion systems.u In those studies the Neel temperature 
is found to decrease as the applied pressure increases and 
superconducting order is found to develop in a limited 
range of applied pressures when the Neel temperature is 
reduced below ~1K. 

It is interesting to see how the mean field critical tem- 
peratures in our model vary with the model parameters 
which, in principle, should depend on externally app- 
plied pressure. Increasing pressure should, presumably, 
make the hybridization V and the conduction band hop- 
ping t increase while keeping U and possibly J almost 
constant .Ell In Figure ^ we plot the mean field critical 
temperatures as function of V while keeping the ratio 
V/t fixed. An interesting feature occurs for the d- wave 
pairing symmetry: if the temperature T* that marks the 
onset of m ^ solutions is higher than the supercon- 
ducting temperature T c then d-wave superconductivity 
and local moments coexist at low temperature; but once 
T c becomes greater than T* the latter drops abruptly and 
we cannot find magnetic solutions at the lowest temper- 
atures. At zero temperature we have a quantum phase 
transition that can be tuned using the external pressure: 



as pressure is reduced (at zero temperature) the ground 
state nature of the system changes from non-magnetic 
but superconducting to magnetic and superconducting 
at the critical value V c /eo ~ 0.38. We find this transition 
to be of first order. Such behaviour is not observed for 
s— or p— wave symmetries. 

Going beyond mean field, and in general terms, the 
local moments due to the /-electrons are progressively 
quenched as the temperature lowers. In dilute systems 
the picture is well understood as due to the Kondo screen- 
ing by the conduction electrons. In dense systems how- 
ever the picture is more involved. At low temperatures 
the local moments are not completely quenched. For 
several t/-based materials like URU2S12, UPt 3 or UPdz 
the remaining moments are quite small of the order of 
0.01 — 0.03/ib but for other systems like UPd2Al3 the 
local moment is quite large of the order of 0.85/is. A 
reason is the Kondo-Nozieres compensation theorem that 
states that in the lattice there are not enough available 
conduction electrons to quench the local moments. 

The Ce-based compounds like CeCu 2 Si2 and 
CeCu§Au show competition between magnetism and su- 
perconductivity if doped and/or under pressure with a 
phase higher in temperature where spin fluctuations are 
significant. The behavior of these systems has been 
interpreted as due to the vicinity to a .quantum criti- 
cal pointO. Two pictures arise howeverca: in the first 
one the Kondo temperature is high (the moments are 
quenched at a finite temperature) and when the system 
approaches the quantum critical point there are no free 
moments (assuming that quenching is complete). Then 
the system has to order due to a Fermi surface instabil- 
ity of the spin density wave type. Another possible sit- 
uation is one in which the moments are not completely 
quenched down to T = and arc free to orient themselves 
leading to magnetism. In the case of CeCu§- x Au x J& 
has been recently found that the second picture holdsE3. 
On the other hand the high value of_the Kondo temper- 
ature for the CeCu2<S«2 compoundel indicates possibly 
that the first scenario should hold. Our results show 
that G?-wave pairing excludes magnetism through a quan- 
tum phase transition. The treatment of a spin-density- 
wave would be mathematically similar and also lead to 
the exclusion of the spin-density-wave by d-wave pair- 
ing as pressure increases. This means that our results 
might be of relevance for understanding the interplay 
between superconductivity and magnetic order in the 
d-wetve heavy- fermion compound CeCu2Si2- This sys- 
tem exhibits a magnetic "A phase" at low temperature 
whose detailed nature is not yet known. Increasing pres- 
sure reduces th£_critical temperature Ta of the A phase. 
Recent studiesalHrtJ of CeCu2Si2 samples near stoichio- 
metric composition have shown that a d-wa.ve supercon- 
ducting phase expels the magnetic "A phase" when Ta 
approaches T c under increasing pressure. Our results in 
Figure || show that d— wave superconductivity destroys 
the local moments when T* meets T c as function of in- 
creasing hybridization. 
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FIG. 2. Mean field magnetic (T*) and superconducting 
(T c ) critical temperatures versus V/|eo|, for a constant ra- 
tio V/t = 0.5. The parameters are U = 2.8, J = -1.7, 
eo = —3.4, and n = 1. The various curves for the magnetic 
transition coincide up to V/eo ~ 0.38. 

In summary, we have studied the formation of local mo- 
ments in the periodic Anderson model in the presence of 
superconducting correlations. We identify the local mo- 
ment regime with the symmetry-breaking saddle-point of 
the effective action for the /-sites in a Hartree-Fock de- 
coupling of local Coulomb repulsion. We have found that 
local moments tend to order antiferromagnetically. We 
note that this does not imply that antiferromagnetic or- 
der arises in a Slater scenario. Magnetic order in heavy 
fermion systems is actually believed to be due to the ef- 
fective RKKY interaction between local moments. This 
means that once we know the region in parameter space 
where local moments form, further studies of magnetic 
order, Kondo screening, etc, possibly require ap-,effective 
interaction Hamiltonian for the local moments£3 By sim- 
ulating the effect of increasing pressure, we have found 
that d-wave superconductivity competes (rather than co- 
exists) with local moment formation at low temperature 
above some critical pressure. 
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